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Administrative control of the Bell System telephone network re- 
quires demand- servicing procedures that maintain network service 
in a cost-effective manner and an administrative measurement plan 
that determines whether the network provisioning process is provid- 
ing a proper balance between network service and network utiliza- 
tion. The traffic measurements available are busy-season, busy-hour, 
trunk-group usage and average blocking. Statistical variation of 
these measurements can lead to an incorrect assessment of network 
performance and an unstable network provisioning process, resulting 
in costly excess reserve capacity in the network. To minimize sensi- 
tivity to volatile measurements, administrative measurement bands 
have been employed to define a new demand- servicing procedure and 
a new trunk service measurement plan for use within the Bell System. 
This paper derives the probability distribution function for measured 
blocking, which is required to construct the measurement bands, and 
outlines the new Bell System administrative procedures. The distri- 
bution has a discrete component, the probability of no blocking, and 
a continuous component. The last function is well approximated by 
a two-parameter beta distribution. The resulting measurement bands 
are a strong function of trunk group size and a weak function of 
traffic characteristics. Small trunk groups (fewer than 10 circuits) are 
much more volatile than large groups. Large groups, properly sized 
for 0.01 blocking, should not experience more than a 0.03 average 
blocking during the 20-day busy season busy hour (bsbh), although 
a properly sized two- circuit group may experience as much as 0.08 
average bsbh blocking. 

I. INTRODUCTION 

The traffic measurements available for network administration are 
also used to estimate the number of trunks required to meet the 
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network engineering objective: an average blocking of 0.01 on the final 
trunk groups and a minimum-cost configuration for the high-usage 
trunk groups. 

The statistical variation of the traffic measurements during the 
20-day busy season causes volatility in the trunks-required estimates. 
This can lead to an incorrect assessment of network performance and 
an unstable network provisioning process. Such problems eventually 
result in unnecessary excess reserve capacity in the trunk and facility 
network, which can exceed 10 percent of that actually required. Ac- 
cordingly, a stable provisioning process is important to minimize the 
reserve capacity while maintaining network service at the objective 
level. 

To minimize the sensitivity of the provisioning process to volatile 
measurements, a concept of administrative measurement bands has 
been used recently to define a new demand-servicing procedure and a 
new trunk service measurement plan for use within the Bell System. 
The demand-servicing procedure provides specific guidelines for de- 
tecting and correcting service problems. The trunk service measure- 
ment plan provides a measure of the balance between network service 
and network utilization. 

The probability distribution functions for the measured average 
blocking, B, and estimated trunks required, c, are required to define 
the administrative measurement bands. The distribution for c is pro- 
vided in Refs. 1 and 2. The purpose of this paper is to establish the 
distribution of B and outline the new Bell System administrative 
procedures. 

The distribution for B has a discrete component, the probability of 
no blocking, and a continuous component. The continuous component 
is well approximated by a two-parameter beta distribution. 

The distribution for B and the previously derived distribution for c 
were used to obtain probability intervals, called measurement bands, 
for both B and c. These bands have been used to establish a new 
demand-servicing policy and a new trunk service measurement plan 
for the Bell System. 

The demand-servicing procedure provides specific guidelines for 
detecting and correcting service problems. For the larger final trunk 
groups, demand servicing is not initiated unless B exceeds 0.03. For 
the smaller groups, the threshold is a function of trunk-group size, and 
increases as the number of circuits decreases. 

The trunk-service measurement plan is based on a comparison of 
the network-sample distribution of B and c, with the corresponding 
theoretical probability distributions. The deviation between the sam- 
ple and theoretical distributions is measured and reduced to a single 
performance index which is a balanced measure of both service and 
utilization. 
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The probability distribution of B is derived in Section II. Section III 
gives numerical results which illustrate the statistical properties of B, 
and Section IV discusses the major application of the results for 
network administration. 



II. THE DISTRIBUTION FUNCTION FOR AVERAGE BLOCKING 

In the Bell System, the trunk-group blocking probability for final 
trunk groups is estimated by measuring the fraction of blocked calls in 
each of the time-consistent busy hours of the 20 consecutive business 
days in the busy season and then averaging the 20 measurements to 
obtain the busy-season, busy-hour, average blocking. The average 
blocking can be highly variable. The variation is a function of the 
length and number of measurement intervals, the traffic characteristics 
(load, peakedness, and day-to-day load variation), and trunk group 
size. In this section, we establish an approximation for the probability 
distribution for average blocking which is adequate for engineering 
purposes. 

2.1 Mathematical model 

We consider a service system with c servers having exponentially 
distributed service times and serving traffic under a blocked-calls- 
cleared service discipline. The mean service time is the time unit. The 
system is observed during n disjoint measurement intervals Ji, • • • , 
./„, each of length T. During Jj, the interarrival times are independent 
and identically distributed (i.i.d.) with mean 1/cy. The peakedness z of 
the offered traffic is constant over all the intervals.* The system is in 
statistical equilibrium during each interval, and the initial point of 
each interval is a stationary (random) point for the arrival process. 
The loads ai, • • ■ , a* are independent and identically distributed 
according to the distribution function T{a\d, Vd) with mean a and 
variance i>d(the day-to-day source-load variance). 

We use Aj(T) and Oj(T), respectively, to denote the number of 
arrivals (call attempts) and the number of calls which were blocked 
(the overflows) during Jj\ the (measured) observed blocking is Bj = 
Oj(T)/Aj(T). The sample average of the observed blocking is 

1 " 

Bn'-ZBj. 

nj-j 

In this section, we derive an approximation for the probability distri- 
bution function for B„. 



* It is known that the peakedness of the traffic overflowing typical high usage trunk 
groups is relatively insensitive to fluctuation in the day-to-day loads offered to those 
groups (see Ref. 4). 
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The distribution can be expressed as a sum of two components; for 
< b < 1, 

P{B n < b) = P{B n = 0} + P{B n > 0}P{B n < b \S„ > 0}. 

In the next subsection, we derive an approximation for the discrete 
component P{B n = 0} and obtain the relations necessary for the 
computation of the conditional distribution. 

2.1.1 The discrete component 

Since a it • • • , a n are i.i.d., B\, • • • , B n are also i.i.d. Thus, 

P{B n = 0} = (P{B l = 0}) n . (1) 

Since B\ = if and only if there are no overflows in J\, we need to 
determine the probability distribution function G c for the length of 
time to the first overflow after the initial point of Ji. After G c is 
determined, we then have (setting ai = a for notational simplicity) 



P{B: = 0}= P{B! = 0|a})dr(a|a, V d ) 



[l-G c (T\a)]dT(a\a, v d ). (2) 



If H c is the conditional distribution for interoverflow times given the 
offered load a, then G c is the corresponding excess-time distribution; 5 
i.e., 



G c (T|a) = a I [I - H c (t\a)]dt, 
Jo 



(3) 



where ao is the mean overflow load corresponding to the offered load 
a. Hence we need to specify the interoverflow distribution H c . 

The Laplace-Stieltjes transform of H c is easily obtained from Ref. 5. 
However, the analytical inversion of the transform is not practical. If 
necessary, one could do the inversion numerically. However, for tele- 
traffic applications, it is more efficient to approximate the interover- 
flow distribution as a mixture of exponentials using the Interrupted 
Poisson Process (ipp) 6 with a three-moment match. The approximation 
is 

H c (t | a) « 1 - kie- rtt - k 2 e~ r2t (4) 

where the parameters k\, k-i, r x and r-i are functions of the load ao and 
the peakedness zo of the overflow traffic. The required relations are 
given in Ref. 6. 
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It now follows from eqs. (3) and (4) that 

1 - G c (T\a) s 1 - (J- (1 - e- 7 ) + - (1 - e^ 7 ) 
n r 2 



From Ref. 6, we know that 



ri r 2 oto 



so 



1 - G C (T | a) * « (^ <r r ' r + £ e"^ . (5) 



Combining eqs. (1), (2), and (5), we have 
P{B n = 0} = 



ab|-c^ 7, + -^ r )dr(a|a f i> < ,) 



n r 2 



(6) 



where it is understood that ao, Ah, n, and r 2 are functions of a, 2, 
and c. 

The distribution function T(a | a, Vd) of the daily loads is assumed to 
be a gamma distribution. 2 Following Ref. 3, we assume that the source- 
load variance v<t is of the form 

v d = 0.13(a)*-—, (7) 

where $ is an adjustable parameter. For Bell System applications, the 
standard values for $ are 1.5 (low day-to-day variation), 1.7 (medium), 
and 1.84 (high). To illustrate the results, we assume throughout this 
study that the mean holding time is 225 seconds and that the mea- 
surement interval is one hour.* Consequently, the relative measure- 
ment time is T = 16. Results for other cases can be generated easily. 
From (6), it follows that the discrete component decreases as T 
increases. More generally, the spread of the distribution of B n will 
decrease as T increases. The integration in (6) can be carried out 
numerically using a 51-point compound Simpson's rule. 

2. 1.2 The continuous component 

Using a simulation, we found that the conditional distribution 
P {B n < b\B n > 0} could be well approximated over the regions of 
engineering interest by a beta distribution, f 7 We call its mean n$ = 
E(B n | B n > 0) and its variance o} = Var(B„ | B n > 0). 



* By numerical experimentation, we found that the results were not very sensitive to 
moderate deviations in mean holding time. A period of 225 seconds is a Bell System 
average holding time. 

f We first tried a lognormal distribution. However, further experimentation indicated 
that a beta distribution would provide a better fit, particularly near the origin. 
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Parameters. It is straightforward to show that 

I,& l fi " >0,= Wo) (8) 

and, using the analog of (8) for higher moments than the first, 

Var(J3„ | B n > 0) = J^f'l - E\B n \ B n > 0)P{B n = 0} . (9) 

Hence, only the ordinary moments of B n need be specified. An ap- 
proximation for the mean was given in Ref. 3; namely, 



E (B n ) = E(Bi | a)aT(a \ a, v d ) (10) 

Jo 

where 

f VaxA(T) Cov[A(T), O(T)]} t „. 

and B(c, a, z) is the trunk group call-congestion corresponding to the 
load a and peakedness z. Algorithms for the computation of the 
(conditional) moments of A(T) and 0(T) (given a and z) are given in 
Ref. 5. 

Since B u • • • , B n are i.i.d., 

Var(B„)=-Var(£i). (12) 

n 

Using the methods described in Ref. 2, Var(Bi) can be expressed as 
Var(Si) = £(Var(fli|a)) + Var^B^a)) 



Var(Bi | a) dT(a \ a, v d ) 

i-ao 

+ I E 2 {B X | a) dT(a \ a, v d ) - ^ 2 (Bi), (13) 

Jo 

where Var(Bi | a) is computed using the algorithms derived in Ref. 5. 
Computation. As noted above and illustrated below, simulated data 
have shown that the conditional distribution P{5„<6|B„>0}is 
well approximated by a beta distribution h{8, tj). That is, for 

0<6<1, 

•b 

x s ~ l (l - x)" -1 dx 

P{B n < b | B n > 0} s I b (8, V ) = ^ , (14) 

x s -\l - X)" -1 dx 

Jo 
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where (see Ref. 7, p. 930) 5 and rj are related to fip and a\ (defined 
above) by 

' . 2 *I 

M/? = P , and a% = — 2 



<5 + 7j p (5 + 7]) 2 (5 + ij + 1) " 

The corresponding inverse relations are 

5=(1-M^-^ (15) 

and 

t>=(1-m((1-/^)(^)-i). (16) 

The integral in the denominator of (14) is the beta function 38(8, tj) 
and is related to the gamma function T(x) as follows (see Ref. 7, page 
258): 

@(8, tj) = x s ~ l (l - j)"" 1 dx (17) 

Jo 

- ™ r <"> (18) 

r<« + ,) • (18) 

To cover the range of engineering interest, it is necessary to use 
logarithms to carry out the computation. Standard programs are 
available for computing the logarithm of the gamma function. 

2.2 Comments 

Collecting the preceding results from Section II, we have 

P{B n < b) = P{B n = 0} + P{B n > 0}P{B„ <b\B„> 0}. 

The approximations for the individual terms are as follows: From 
Section 2.1.1, 



P{B n = 0) = 



\( a J^ e -nT + ^ e -r 2 T\ dr{al - Vd) 



where a , k\, k 2 , n and r> are functions of a, z, and c as given in Ref. 
6. From Section 2.1.2, 

P{B n < b\B n >0) ~ h(8,-n) 

n« + ti) C . , 

1 x s ~ l (l - x)"' 1 dx, 



r(8)T(r,) J 
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where 8 and tj are related to the conditional mean and variance of B n 
as given in eqs. (15) and (16). The integrations are carried out using a 
51-point compound Simpson's rule. 

2.3 Quantiles 

The p-quantile of a continuous distribution F is the value £ p such 
that F{£ p ) = p; i.e., £ P = F~ 1 (p), where F _1 denotes the inverse function 
forF. 

From the preceding section, we have 

P{Bn < b) » P{B n = 0} + P{Bn > 0}I b {8, -n). 

Using this result, the quantiles of the distribution are easily expressed 
in terms of the quantiles of the beta distribution. That is, for a given 
p between P{B n = 0} and 1, the value £ p such that 

is given approximately by 

£p * I&» (8, t,) (19) 

where 

p-P{B n = 0) 



q(p) = 



P{B n >0} 



and Iq l (8, rj) denotes the g-quantile of the beta distribution (of which 
I' 1 is the inverse). Standard programs are available for computing the 
beta distribution and its quantiles. 8 

III. NUMERICAL RESULTS 

In this section we present a few examples of the results which have 
been obtained from the approximations developed in Section II. We 
consider the discrete component and the cumulative distribution of 
B n , the corresponding probability density function, and a set of prob- 
ability intervals. Applications of the results are described in Section 
IV. 

3. 1 The discrete component 

The discrete component, P{B n = 0}, has several potential engineer- 
ing applications in addition to being a fundamental part of the distri- 
bution of B n . Accordingly, in this subsection we verify the accuracy of 
the approximation for P{B n = 0} and illustrate the general character- 
istics which should be useful in future applications. 

Approximation (6) was computed for a wide range of system param- 
eters and compared with estimates obtained from a random-walk 
simulation of a trunk group satisfying the assumptions given in Section 
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2.1. Typical results are shown in Fig. 1, where a graph of P{B\ = 0} 
versus trunk group size c is shown for trunk groups engineered at a 
mean blocking of 0.01. The three curves correspond to z = 1 and low 
day-to-day variation (<f> = 1.5), 2 = 2 and medium variation (</> = 1.7), 
and 2 = 4 and high variation (<£ = 1.84). 

The agreement between the simulation and the approximation is 
quite good for z = 1 and 2. For 2 = 4, the approximation is a little low, 
but even in this case (which is relatively extreme for most traffic 
applications), the relative difference is less than 7 percent and should 
be adequate for engineering applications. (The approximation was not 
tested for z > 4.) Similar accuracy was obtained in the other cases 
tested, particularly when P{Bi = 0} was on the order of 0.1 or larger 
(i.e., low-blocking applications). The accuracy is statistically more 
difficult to check for smaller values. However, when P{Bi = 0} is 
much less than 0.1, it has very little influence on the distribution of 
B„. Also, notice that, as n increases, P{B n = 0} decreases rapidly; i.e., 
it has very little engineering significance for n larger than 10. 

Figure 2 shows P{B\ = 0} versus the load a for c = 40 trunks and 
several values of peakedness (no day-to-day variation). As the figure 
shows, P{B\ =0} is a monotone decreasing function of both a and z 
when no day-to-day variation is present. Similar results were also 
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Fig. 1 — P{B\ = 0} vs trunk-group size. 
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Fig. 2— P{Bi = 0} vs load for C = 40 trunks. 

obtained for any fixed level of variation. Note that the separation of 
the curves introduces the possibility of estimating a and z by using 
measurements of P{Bi = 0} and usage. 

The effect of day-to-day variation for the case z = 1 is also illustrated 
in Fig. 2. Day-to-day load variation decreases P{B X = 0} at lower loads 
and increases it at higher loads. At a high mean load, day-to-day 
variation increases the likelihood of a low load and therefore increases 
P{B\ = 0}. At a low mean load, day-to-day variation increases the 
probability of a relatively larger load and thereby decreases 
P{B\ = 0}. The difference can be considerable, and day-to-day varia- 
tion should be taken into account when computing the discrete com- 
ponent. 

3.2 The cumulative distribution 

The approximation for the cumulative distribution 

P{B„ < b) « P{B n = 0} + P{B n > 0} h(S, t,), < b < 1, 

(derived in Section II) was computed for several system parameters 
and compared with estimates of the distribution obtained from the 
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simulation. Typical results are shown in Fig. 3 for c = 40 trunks, 
z — 1, <}> — 1.5 (low variation), and n = 1, 5, and 20 days. The 
approximation compares quite well with the simulated data. 

A number of additional cases were generated to determine the range 
of validity of the approximation. We considered c = 3, 5, 10, 40, and 
100, 2=1 with low day-to-day variation and 2 = 4 with high variation. 
For all combinations of load and peakedness, the accuracy of the 
approximation increases with c and should be adequate for engineering 
applications with c > 10. For the smaller values of c, c < 5, the 
approximation did not agree as well with the simulated data. In those 
cases, a simulation is probably the most appropriate tool to estimate 
the distribution. 



3.3 The probability density 

Since the distribution of B n has a discontinuity at zero, the proba- 
bility density function is not defined at zero in the usual sense; i.e., the 
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Fig. 3 — Distribution of observed blocking. 
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density function contains an impulse at b = 0. 9 For b > 0, an approxi- 
mation for the probability density is obtained by differentiating the 
approximate distribution to obtain 

/■(fe) = iL P{B n < 6} = P{B„ >0} 



rf& 



#(«, i?) 



where 5, tj and ^ were defined in (15) to (17). 

The number of hours of data has a strong influence on the shape of 
the density function. Figure 4 illustrates the cases n = 1, 5, and 20. For 
each case, c = 30, z = 1, the day-to-day variation is low, and the 
average blocking is 0.01. Notice that the mode of each density is 
located to the left of the mean value 0.01. This is probably caused by 
the nonlinear effect on trunk-group blocking of a variable load. 

IV. ADMINISTRATIVE MEASUREMENT BANDS 

The statistical variation of the traffic measurements during the 20- 
day busy season causes volatility in the measured average blocking 
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Fig. 4 — Probability density for observed blocking. 
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and in the trunks-required estimates. Such volatility can cause unnec- 
essary demand servicing which will result in excess reserve trunk 
capacity. Unless the provisioning process is stabilized, this excess 
capacity can eventually exceed 10 percent of the actual required 
trunks. 

A concept of administrative measurement bands has been used to 
define a stable demand servicing procedure and trunk administrative 
measurement plan. In this section, we describe the measurement bands 
and summarize their planned use in the Bell System provisioning 
process. 

4. 1 Measured blocking 

The numerical results presented in Section III show that when a 
trunk group is properly engineered for an average blocking of 0.01, the 
observed average blocking may deviate considerably from 0.01. In fact, 
the probability density function for B n is skewed to the right suffi- 
ciently that the observed blocking will fall below the 0.01 average 
about 60 percent of the time. Consequently, if a naive observer took a 
single look at a properly engineered network, he might conclude that 
many of the probability-engineered trunk groups were overprovided 
(observed average blocking somewhat less than 0.01). On the other 
hand, nearly 40 percent of the groups would experience blocking in 
excess of 0.01 and thereby be of concern to the trunk engineer who 
sized the groups for an 0.01 objective. 

To account properly for the volatility of measured blocking, intervals 
of acceptable and anticipated measured average blocking have been 
calculated for Bell System network administration. The resulting 
intervals are illustrated in Fig. 5. 

The analytical approximation was used to generate the curves for 
c > 7. A simulation was used to provide the results for 2 < 
c < 6. For c < 5, the analytical approximation is not adequate for our 
purposes. Moreover, a simulation is inexpensive and effective for the 
small trunk-group sizes. 

The upper and lower curves correspond to the theoretical 0.99- and 
0.01-quantiles, respectively. However, simulation results have shown 
that these theoretical quantiles are the effective end points of the 
distribution. The lower curve vanishes for c < 50 because the smaller 
groups can be properly sized and still, with probability greater than 
0.01, experience no busy-season, busy-hour blocking. The two central 
curves bound an 80-percent probability interval. That is, 80 percent of 
all the average-blocking measurements should lie within the interval 
if the network is properly sized. 

The curves are fairly constant for c > 10. However, they increase 
substantially for the smaller groups, illustrating the statistical volatility 
of traffic measurements on the smaller groups. 
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Fig. 5— Trunk service measurement plan, average blocking measurement bands. 



4.2 Trunks required 

The distribution function for estimates of trunks-required is pro- 
vided in Refs. 1 and 2. It has been used to construct administrative 
bands that are analogous to those for average blocking. The results for 
low day-to-day variation are given in Fig. 6. 

The vertical axis represents acceptable deviations in the estimates. 
For example, if a 10-trunk group were properly sized, then the esti- 
mates should fall between 9 and 11 about 80 percent of the time. The 
results shown must be modified for medium and high day-to-day 
variation; the bands increase in width as the level of variation increases. 

4.3 Network administration 
4.3. 1 Demand servicing 

Trunk network administration is composed of two major operations, 
demand servicing and planned servicing. Planned servicing uses the 
forecasts of demand in the next busy season to decide how many and 
where circuits are required to insure adequate network service. De- 
mand servicing refers to the use of trunk group measurements to 
detect and correct current service problems. The measurement bands 
for & and c have been used to establish a new demand servicing policy 
which is planned to be used in the Bell System. 10 

Based on the new policy, a network cluster is assumed to be 
providing adequate service provided the average blocking for the final 
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trunk group does not exceed the upper curve. If the blocking exceeds 
the threshold, then the trunks-required administrative bands are used 
to determine cost-effective trunk additions which will restore network 
service. 

4.3.2 Measurement plan 

A new Trunk Administrative Measurement Plan was recently issued 
for use within the Bell System. The trunking accuracy component of 
the plan, called the Trunk Service Measurement Plan (tsmp), uses the 
measurement bands for B and c to measure the balance between 
network service and network utilization. 

The tsmp compares the network-sample distribution of B and c with 
the theoretical distribution; a 5-point comparison is made using the 0- 
10-80-10-0 intervals illustrated in Figs. 5 and 6. The deviation between 
the sample and theoretical distributions is measured and reduced to a 
single performance index which is a measure of the balance between 
network service and utilization. 

V. ACKNOWLEDGMENTS 

P. J. Burke (Bell Laboratories) used the measurement band concept 
to establish the Trunk Administrative Measurement Plan. C. R. Szelag 
(Bell Laboratories) constructed the new Trunk Demand Servicing 
Policy. 

BLOCKING DISTRIBUTIONS 843 



REFERENCES 

1 S R. Neal and A. Kuczura, "A Theory of Traffic-Measurement Errors for Loss 

Systems With Renewal Input," B.S.T.J., 52, No. 6 (July- August 1973), pp. 967- 

990. 

2 D W. Hill, "The Effects of Load Variation on Measurement Accuracy," B.S.T.J., 

56, No. 4 (April 1977), pp. 561-573. 

3 D W. Hill and S. R. Neal, "Traffic Capacity of a Probability-Engineered Trunk 

Group," B.S.T.J., 55, No. 7 (September 1976), pp. 831-842. 

4. M. Berger and S. R. Neal, "A Sensitivity Study of Traffic-Parameter Estimation 

Procedures Used for Engineering Trunk Groups," Proceedings of the Seventh 
International Teletraffic Congress, Stockholm, June 1, 1973, paper 542. 

5. A. Kuczura and S. R. Neal, "The Accuracy of Call-Congestion Measurements for 

Loss Systems with Renewal Input," B.S.T.J., 51, No. 10 (December 1972), pp. 
2197-2208. 
6 A Kuczura, "The Interrupted Poisson Process As An Overflow Process, B.S.T.J., 
52, No. 3 (March 1973), pp. 437-448. 

7. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Applied 

Mathematics Series 55, Washington, D.C.: National Bureau of Standards, June 

1964. 

8. C. Hastings, Approximations for Digital Computers, Princeton, N.J.: Princeton U. 

Press, 1955, pp. 191-192. 

9. M. J. Lighthill, An Introduction to Fourier Analysis and Generalized Functions, 

Cambridge, England: Cambridge U. Press, 1959, pp. 16-17. 
10. C. R. Szelag, "Trunk Demand Servicing in the Presence of Measurement Uncer- 
tainty," B.S.T.J., this issue, pp. 845-860. 



844 , THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1 980 



